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ABSTRACT 

From Space Telescope Imaging Spectrograph (STIS) observations of Uranus in 2012 , when good 
views of its north polar regions were available, we found that the methane volume mixing ratio 
declined from about 4% at low latitudes to about 2% at 60° N and beyond. This depletion in the 
north polar region of Uranus in 2012 is similar in magnitude and depth to that found in the south polar 
regions in 2002. This similarity is remarkable because of the strikingly different appearance of clouds 
in the two polar regions: we have never seen any obvious signs of convective activity in the south polar 
region, while the north has been peppered with numerous small clouds thought to be of convective 
origin. Keck and Hubble Space Telescope imaging observations close to equinox at wavelengths of 1080 
nm and 1290 nm, with different sensitivities to methane and hydrogen absorption but similar vertical 
contribution functions, imply that the depletions were simultaneously present in 2007, and at least 
their gross character is probably a persistent feature of the Uranus atmosphere. The depletion appears 
to be mainly restricted to the upper troposphere, with the depth increasing poleward from about 30° 

N, reaching ^4 bars at 45° N and perhaps much deeper at 70° N, where it is not well constrained 
by our observations. The latitudinal variations in degree and depth of the depletions are important 
constraints on models of meridional circulation. Our observations are qualitatively consistent with 
previously suggested circulation cells in which rising methane-rich gas at low latitudes is dried out 
by condensation and sedimentation of methane ice particles as the gas ascends to altitudes above the 
methane condensation level, then is transported to high latitudes, where it descends and brings down 
methane depleted gas. Since this cell would seem to inhibit formation of condensation clouds in regions 
where clouds are actually inferred from spectral modeling, it suggests that sparse localized convective 
events may be important in cloud formation. A more complex meridional circulation pattern may 
be necessary to reproduce the observed cloud distribution, but microwave observations appear to be 
most compatible with a single deep circulation cell. The small-scale latitudinal variations we found in 
the effective methane mixing ratio between 55° N and 82° N have significant inverse correlations with 
zonal mean latitudinal variations in cloud reflectivity in near-IR Keck images taken before and after 
the HST observations. If the CH 4 /H 2 absorption ratio variations are interpreted as local variations 
in para fraction instead of methane mixing ratio, we find that downwelling correlates with reduced 
cloud reflectivity. While there has been no significant secular change in the brightness of Uranus at 
continuum wavelengths between 2002 and 2012 , there have been significant changes at wavelengths 
sensing methane and/or hydrogen absorption, with the southern hemisphere darkening considerably 
between 2002 and 2012 , by ^25% at mid latitudes near 827 nm, for example, while the northern 
hemisphere has brightened by ^25% at mid latitudes at the same wavelength. 

Subject headings: Uranus, Uranus Atmosphere; Atmospheres, composition; Atmospheres, dynamics 


1. INTRODUCTION 

Uranus experiences the solar system’s largest fractional 
seasonal forcing because its spin axis has a 98° inclina¬ 
tion to its orbital plane. It is thus not surprising to see a 
time-dependent north-south asymmetry in Uranus’ cloud 
structure. An especially interesting asymmetry noted by 
iSromovskv et al.l (| 2012 [ ) is the continued complete ab- 
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sence of discrete cloud features south of 45° S, while nu¬ 
merous discrete cloud features have been observed north 
of 45° N in recent near-IR H-filter Keck images (Fig. [1]). 
Voyager imaging in 1986 recorded no bright cloud fea¬ 
tures between the south pole and 45° S. Nor were any 
seen in near-IR Hubble Space Telescope (HST) images 
taken with the Near-IR Camera and Multi-Object Spec¬ 
trometer (NICMOS) beginning in 1997, nor in Keck im¬ 
ages beginning in 2003. Some mechanism appears to be 
inhibiting convection at high southern latitudes that is 
not present at high northern latitudes. A very significant 
and possibly related result came from analysis of 2002 
Space T elescope Imaging Spectrometer ( STIS) observa¬ 
tions by iKarkoschka and Tomaskol (|2009[) , subsequently 
re ferenced as KT2 0 09, a nd confirmed by the analysis 
of ISromovskv et ^ (| 2011 [) . They found a strong deple¬ 
tion of methane in the upper troposphere (down to a few 
bars at least) at high southern latitudes, suggesting a 
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downwelling flow at these latitudes, which would tend to 
inhibit convective cloud formation. This raised the pos¬ 
sibility of a connection between methane depletion and 
a lack of discrete cloud features, suggesting that high 
northern latitudes, where discrete clouds are seen, might 
not be depleted in methane. If so, the methane deple¬ 
tion might be a seasonal effect. In 2002, five years before 
equinox, the sub-observer latitude was 21.4° S and the 
north polar region was not visible, so that testing this hy¬ 
pothesis would require new observations when the north 
polar region was exposed to view. We thus proposed new 
STIS observations of Uranus, which were obtained in late 
September 2012 (HST program 12894, Sromovsky, PI), 
five years after equinox, when the sub-observer latitude 
was 19.5° N. The analysis of these new observations is 
the primary topic in what follows. 

Constraining the mixing ratio of CH 4 on Uranus is 
based on differences in the spectral absorption of CH 4 
and H 2 , illustrated by the penetration depth plot of Fig. 
EK. Methane absorption dominates at most wavelengths, 
but hydrogen’s Collision Induced Absorption (CIA) is 
relatively more important in a narrow spectral range near 
825 nm. Model calculations that don’t have the correct 
ratio of methane to hydrogen lead to a relative reflectivity 
mismatch near this wavelength. An example is shown in 
Fig. EB-D, in which model calculations are compared to 
2002 STIS observations assuming methane profiles with 
2.2% and 4.0% deep volume mixing ratios. The two as¬ 
sumptions lead to very different errors in the vicinity of 
825 nm, clearly i ndicating that the larger mixing ratio 
is a better choice. iKarkoschka and Tomaskol ()2009f ) used 
this spectral constraint to infer a methane mixing ratio 
of 3.2% at low lati tudes, but dropping to 1.4% at high 
southern latitudes. iSromovskv et al.l (j 201 l[ ) analyzed the 
same data set, but used only temperature and mixing 
ratio p rofiles that were consistent with the iLindal et al.l 
(jl987[ ) refractivity profiles. They confirmed the depletion 
but inferred a somewhat higher mixing ratio of 4% at low 
latitudes and found that better fits were obtained if the 
high latitude (down to ^2-4 bars). Subsequently, 2009 
groundbased spectral observations at the NASA Infrared 
Telescope Facility (IRTF) using the SpeX instrument, 
which provided cov erage of th e key 825-nm spectral re¬ 
gion, were used by iTice et M ] (imi to infer that both 
polar regions were weakly depleted in methane, but they 
inferred lower methane mixing ratios, smaller latitudinal 
variations, and higher u ncertainties than the STI S-based 
analysis of KT2009 and ISromovskv etld] (j 201 l[ ). These 
lower IRTF-based values might be a result of lower spa¬ 
tial resolution combined with worse view angles into the 
polar regions than obtained by HST observations. 

In the following we begin with a description of our 
2012 HST/STIS observations, and the complex data re¬ 
duction and calibration procedures. We then describe 
the relatively direct implications of latitudinal spectral 
variations and simplified model results. Finally, we de¬ 
scribe our radiation transfer calculation methods and the 
results of constraining the cloud structure and methane 
distributions to fit the observed spectra. We show that 
the methane depletion is indeed present in both polar 
regions but of significantly different character. The re¬ 
sults provide no obvious explanation for the asymmetry 
in polar cloud structure, but do raise some important 
questions about cloud formation. The important impli¬ 


cations for meridional circulation models are discussed 
before summarizing our results. 

2. OBSERVATIONS 

Our 2012 observations used four HST orbits, three of 
them devoted to STIS spatial mosaics and one orbit to 
Wide Field Camera 3 (WFC3) support imaging. The 
STIS observations were taken on 27-28 September 2012 
and the WFC3 observations on 30 September 2012 . Ob¬ 
serving conditions and exposures are summarized in Ta¬ 
ble [H 

2 . 1 . STIS spatial mosaics. 

STIS observations used the G430L and G750L grat¬ 
ings and the CCD detector, which has ^0.05 arcsecond 
square pixels covering a nominal 52"x52" square field 
of vi ew (FOV) and a spectr al range from ^200 to 1030 
nm (jHernandez et al.l 120121 ). Using the 52"x0.1" slit, 
the resolving power varies from 500 to 1000 over each 
wavelength range due to fixed wavelength dispersion of 
the gratings. Observations had to be carried out within 
a few days of Uranus opposition (29 September 2012 ) 
when the telescope roll angle could be set to orient the 
STIS slit parallel to the spin axis of Uranus. 

One STIS orbit produced a mosaic of half of Uranus us¬ 
ing the CCD detector, the G430L grating, and 52" x0.1" 
slit. The G430L grating covers 290 to 570 nm with a 
0.273 nm/pixel dispersion. The slit was aligned with 
Uranus’ rotational axis, and stepped from the evening 
limb to the central meridian in 0.152 arcsecond incre¬ 
ments (because the planet has no high spatial resolution 
center-to-limb features at these wavelengths we used in¬ 
terpolation to fill in missing columns of the mosaic). Two 
additional STIS orbits were used to mosaic the planet 
with the G750L grating and 52" x0.1" slit (524-1027 nm 
coverage with 0.492 nm/pixel dispersion), with limb to 
central meridian stepping at 0.0569 arcsecond intervals. 
This was the same procedure that was used successfully 
for HST program 9035 in 2002 (E. Karkoschka, P.I.). As 
Uranus’ equatorial radius was 1.85 arcseconds when ob¬ 
servations were performed, stepping from one step off the 
limb to the central meridian required 13 positions for the 
G430L grating (at 0.152 arcseconds/step) and 36 for the 
G750L grating (at 0.0569 arcseconds/step). Two orbits 
were needed to complete the G750L grating observations, 
spanning a total time of 2 h 17 m, during which Uranus 
rotated 47°. This rotation was not a problem because of 
the high degree of zonal symmetry of Uranus and because 
our analysis rejected any small scale deviations from it, 
such as rare discrete cloud features. 

Exposure times were similar to those used in the 2002 
program, with 70-second exposures for G430L and 86 - 
second exposures for G750L gratings, using the 1 elec- 
tron/DN gain setting. These exposures yielded single¬ 
pixel signal-to-noise ratios of around 10:1 at 300 nm, 
> 40:1 from around 400 to 700 nm, and decreasing to 
around 20:1 (methane windows) to < 10:1 (methane ab¬ 
sorption bands) at 1000 nm. 

2 . 2 . Supporting WFC3 imaging. 

Since STIS images can be radiometrically calibrated 
for point sources or infinitely-extended sources, and 
Uranus is neither, an empirically determined correction 
function must be applied to the images as a function of 
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Fig. 1.— Shift-and-add high-pass filtered projections of Uranus’ 2003 South (L) and 2011 North (R) polar regions reveal a large asymmetr y 
in cloud features poleward of 45°, suggesting very different degrees of polar convection, a likely seasonal effect (ISromovskv et ani2012|j . 
The images were obtained with the Keck 2 telescope and N1RC2 instrument using an H filter (1.62 fim central wavelength). The small 
bright features in the polar region in the 2011 image are 1-2% brighter than their surroundings. The features in the lower right quadrant 
of that image are heavily distorted by projection artificats. 


TABLE 1 

Science exposures from HST program 12894. All STIS spectra used the 52" x 0.1" slit. 


Orbit 

Start 

Date (UTH) 

Start 

Time (UTH) 

Instrument 

Filter or 
Grating 

Exposure 

(sec) 

No. of 
Exp. 

Phase 
Angle (°) 

1 

2012-09-27 

21:22:29 

STIS 

MIRVIS 

10.1 

1 

0.09 

I 

2012-09-27 

21:38:11 

STIS 

G430L 

70.0 

13 

0.09 

2 

2012-09-27 

22:56:43 

STIS 

G750L 

86.0 

18 

0.08 

3 

2012-09-28 

00:32:26 

STIS 

G750L 

86.0 

18 

0.08 

4 

2012-09-30 

22:44:50 

WFC3 

F336W 

30.0 

1 

0.09 

4 

2012-09-30 

22:46:35 

WFC3 

F467M 

16.0 

1 

0.09 

4 

2012-09-30 

22:48:15 

WFC3 

F547M 

6.0 

1 

0.09 

4 

2012-09-30 

22:49:39 

WFC3 

F631N 

65.0 

1 

0.09 

4 

2012-09-30 

22:52:08 

WFC3 

F665N 

52.0 

1 

0.09 

4 

2012-09-30 

22:54:15 

WFC3 

F763M 

26.0 

1 

0.09 

4 

2012-09-30 

22:56:02 

WFC3 

F845M 

35.0 

1 

0.09 

4 

2012-09-30 

22:57:56 

WFC3 

F953N 

250.0 

1 

0.09 

4 

2012-09-30 

23:04:27 

WFC3 

FQ889N 

450.0 

1 

0.09 

4 

2012-09-30 

23:16:02 

WFC3 

FQ937N 

160.0 

1 

0.09 

4 

2012-09-30 

23:23:05 

WFC3 

FQ727N 

240.0 

1 

0.09 


On September 28 the sub-observer planetographic latitude was 19.5° S, the observer range was 19.0613 AU 
(2.851530x10^ km), and the equatorial angular diameter of Uranus was 3.6976 arcseconds. 


wavelength. This function was determined for the 2002 
STIS observations using WFPC2 images of Uranus taken 
around the same time as the 2002 STIS Uranus obser¬ 
vations. To ensure that this function was determined as 
well as possible for the Cycle 20 observations in 2012, 
and to cross check the extensive spatial and spectral cor¬ 
rections that are required for STIS observations, we used 
one additional orbit of WFC3 imaging at a pixel scale of 
0.04 arcseconds with eleven different filters spread over 
the 300-1000 nm range of the STIS spectra. These WFC3 
images are displayed in Fig. [3l along with synthetic im¬ 
ages with the same spectral weighting constructed from 
STIS spatially resolved spectra, as described in the fol¬ 
lowing section. The filters and exposures are provided in 
Table [TJ 

3. DATA REDUCTION AND CALIBRATION. 


The STIS pipeline processing used at STScI does not 
yield suitably calibrated two-dimensional spectral images 
for an object like Uranus. Considerable additional effort 
was required to reach a final calibration of these data, 
using techniques developed by KT2009 and closely fol¬ 
lowed in the calibration of the 2012 STIS observations. 
Flat-fielded science images, fringe flats, and WAVECALS 
from the STScI STIS data processing pipeline are the in¬ 
puts to a rather extensive post-processing suite, each step 
of which is described in Section 1 of the on-line analy¬ 
sis supplemental file. (WAVECALS are exposures of arc 
sources used for wavelength calibration.) The output is 
a cube of geometrically and radiometrically calibrated 
monochromatic images of Uranus. However, this cube 
needs a final correction based on comparisons with the 
WEC3 support images. 

Each WFC3 image was deconvolved with an appro¬ 
priate Point-Spread Eunction (PSE) obtained from the 
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Fig. 2.— A: Penetration depth vs. wavelength as limited by 
different opacity sources assumin g th e FI methane profile and ab¬ 
sorption models discussed in Sec. 16.II Note the greater importance 
of H 2 CIA near 825 nm. The dark solid curve includes all ab¬ 
sorbers, while the other curves are for sing le absorbers in isolation . 
B: Example model fit for the FI profile of ISromovskv et al.l (1201111 
with a deep methane volume mixing ratio of 4%, and fit errors for 
deep methane mixing ratios of 4% (C) and 2 . 2 % (D). The fit errors 
near 825 nm are a sensitive indicator of the mixing ratio of CH 4 
to H 2 . 

Tiny Tim code of iKristI (|1995[) . To match the spatial 
resolution of the STIS images, the WFC3 images were 
then reconvolved with an approximation of the PSF given 
in the analysis supplemental file. They were then con- 
verted to I/F using header PHOTFLAM values and the 
iColina et al.l (| 19961 ) solar flux spectrum, averaged over 
the WFC3 filter bandpasses. (The PHOTFLAM value is 
used to convert instrument DN/second to flux units.) To 
obtain a disk-averaged I/F, the planet’s light was inte¬ 
grated out to 1.1 equatorial radii and averaged over the 
planet’s cross section in p ixels, computed using NAIF 
ephemerides (|Actonl 119961) and SPICELIB limb ellipse 
model (SPICELIB is NAIF toolkit software used in gen¬ 
erating navigation and ancillary instrument information 
files.) The disk-averaged I/F was also computed for each 
of the STIS monochromatic images, and the filter- and 
solar flux-weighted I/F was computed for each of the 
WFC3 filter passbands that we used. 

By comparing the uncalibrated STIS synthetic disk- 
averaged I/F to the corresponding WFC3 values, we 


constructed a correction function to radiometrically cal¬ 
ibrate the STIS cube. Figure |4] shows the ratios of STIS 
to WFC3 disk-integrated brightnesses, and the quadratic 
function that we fit to these ratios as a function of wave¬ 
length. We weighted the filters by the square root of 
their band widths, and computed an effective wavelength 
weighted according to the product of the solar spectrum 
and the I/F spectrum of Uranus. The RMS deviation 
of individual filters from the calibration curve given in 
Fig. |4]is 3.7%, but only 0.94% when the FQ937N filter 
is omitted from the RMS calculation (but still included 
in the fit). While the FQ937N ratio deviates by 10% 
from the final calibration curve, it is not likely due to a 
problem with the filter calibration, as a comparison with 
the uranian satellite Ariel found it to be consistent with 
the other filters. More likely is that the darker parts of 
the STIS spectrum require slightly different calibrations 
than the brighter parts due to residual errors in Charge 
Transfer Efficiency (CTE) corrections (described in Sec¬ 
tion 1 of the analysis supplement file). Figure |4] also 
shows the function that was used to radiometrically cal¬ 
ibrate the 2002 STIS data by KT2009. We do not know 
whether the difference in functions is due to a true STIS 
instrument change, to spatial variations across the CCD 
(2012 and 2002 observations were acquired at different 
locations on the CCD), or to a difference in filters used 
in the two calibrations (the 2002 data were calibrated 
using images taken by the Wide Field Planetary Camera 
2 (WFPC2) instrument, the predecessor to WFC3). 

As a sanity check on the STIS processing we com¬ 
pared line scans across WFC3 images to the correspond¬ 
ing scans across synthetic WFC3 images created from 
our calibrated STIS data cubes. These generally were 
in excellent agreement for fi values from 0.3 to 0.8, as 
demonstrated by plots given in Section 2 of the on-line 
analysis supplement file. The most significant discrep¬ 
ancy is in the overall I/F level computed for the FQ937N 
filter, a consequence of our calibration curve being 10% 
high for that filter. This can also be seen in the ratio 
images plotted in Fig. [3l 

3.1. Center-to-limb fitting 

Center-to-limb profiles provide important constraints 
on the vertical distribution of cloud particles as well as 
the vertical variation of methane compared to hydrogen. 
Because Uranus has a high degree of zonal uniformity, 
these profiles are fairly smooth functions that can be 
characterized by a small number of parameters, making it 
possible to constrain the profiles accurately, reducing the 
effects of noise and skipping over any small discrete cloud 
features. Because the observations were taken very close 
to zero phase angle these functions are almost perfectly 
symmetric about the central meridian, so that they de¬ 
pend on only one cosine value (observer and solar zenith 
angles are essentially equal). This is also one reason we 
were able to characterize these profiles by measuring only 
one half of the Uranus disk. These fits also facilitate the 
separation of latitudinal variations from those associated 
with view angle variations. 

Before fitting the center-to-limb (CTL) profile for 
each wavelength, the spectral data are smoothed in the 
wavenumber domain to a resolution of 36 cm“^, which 
equals the resolution we use in computing the Raman 
spectrum, but is four times finer than the sampling we 
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Fig. 3.— WFC3 images of Uranus taken on 30 September 2012 (A-F and M-Q) compared to synthetic band-pass filter images (G-L and 
R-V) created from weighted averages of STIS spectral data cubes using WFC3 throughput and solar spectral weighting. The north pole is 
at at the right. Portions of the synthetic images east of the central meridian are obtained by reflection of the images west of the central 
meridian. That is why the single cloud feature appears twice in images where cloud contrast is sufficient to make the cloud visible. The 
STIS data were acquired three days before the WFC3 data, which were taken when the discrete cloud was not visible. The ratio images 
are stretched to make 0.8 black and 1.2 white. A direct comparison of latitude scans at fixed view angles is provided in the on-line analysis 
supplement file. 


use in constraining cloud models. The smoothing is help¬ 
ful in reducing noise at longer wavelengths without de¬ 
grading important spectral features. For each 1° of lat¬ 
itude from 54° S to 85° N, all image samples within 1° 
of the selected latitude and with /i > 0.175 are collected 
and fit to the empirical function 

/(^) = a + 5/i + c//i, (1) 


S to 72.6° N for 2012 observations and 74.5° S to 31.7° 
N for 2002. Unless otherwise noted all our results are 
derived without extrapolation. 

The CTL fits can also be used to create zonally 
smoothed images by replacing the observed I/F for each 
pixel by the fitted value. Results of that procedure are 
displayed in a later section, with further examples shown 
in the analysis supplement. 


assuming all samples were collected at the desired lati¬ 
tude and using the /i value for the center of each pixel of 
the image samples. Sample fits are provided in Fig.[5]and 
in Fig. 2 of the on-line supplement. Most of the scatter 
about the fitted profiles is due to noise, which is often am¬ 
plified by the deconvolution process. Because the range 
of observed /i values decreases away from the equator at 
high southern and northern latitudes, we chose a moder¬ 
ate value of /i=0.6 as the maximum view-angle cosine to 
provide a reasonably large unextrapolated range of 33.6° 


4. DIRECT COMPARISON OF METHANE AND 
HYDROGEN ABSORPTIONS VS. LATITUDE. 

If methane and hydrogen absorptions had the same de¬ 
pendence on pressure, then it would be simple to estimate 
the latitudinal variation in their relative abundances by 
looking at the relative variation in I/F values with lati¬ 
tude for wavelengths that produce similar absorption at 
some reference latitude. Although this idea is compro¬ 
mised by different vertical variations in absorption, which 
means that latitudinal variation in the vertical distribu- 
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Fig. 4.— A: Disk-integrated I/F spectra before final calibration 
(lower curve) and after final calibration (upper curve), where hor¬ 
izontal red bars indicate synthetic I/F values computed from pre¬ 
calibration spectra and horizontal green bars indicate I/F values 
obtained from WFC3 imaging. B: Synthetic band-pass filter disk- 
integrated I/F values (pre-calibration) divided by corresponding 
I/F values obtained from WFC3 measurements (horizontal bars). 
The red (solid) curve and legend coefficients defining it, were ob¬ 
tained by fitting all ratio values including that obtained from the 
FQ937N filter. The black (dashed) curve is the calibration func¬ 
tion obtained from the 2002 data set. C: Ratio of synthetic disk- 
integrated I/F values obtained from calibrated STIS spectra to the 
corresponding WFC3 values. 


imbcorr 

CTL curves at glot 30.00° 


n.R - • 552.0 nm, X 1.0 



0.2 0.4 0.6 0.8 1.0 

PS file: minn_ctLdiag.ps cos(Emission Angle) Apri 4-1 54855-201 4 


Fig. 5. — Sample center-to-limb fits for 30° N, as described in the 
main text. STIS I/F samples and fit lines with uncertainty bands 
are shown for five different wavelengths indicated in the legends. 
The latitude band sampled for these fits is darkened in the inset 
image of the half-disk of Uranus. 

tion of aerosols can also play a role, it is nevertheless 
useful in a semi-quantitative sense. Thus we explore sev¬ 
eral cases below. 

4.1. Direct comparison of key near-IR wavelength scans 

in 2007 

Our first example is based on a comparison of HST 
NICMOS measurements using an F108N filter (centered 


Level at which 2-way OD=1.00 



Fig. 6. — Near-IR penetration depths vs. wavelength com¬ 
pared to filter transmission for FI08N NICMOS and PaBeta 
Keck/NIRC2 filters, which sense similar atmospheric levels in a 
clear atmosphere, but are dominated by different gas absorptions 
(H 2 and CH 4 respectively). 

at 1080 nm), which is dominated by H 2 CIA, with 
KeckII/NIRC2 measurements using a PaBeta filter (cen¬ 
tered at 1290 nm), which is completely dominated by 
methane absorption. [The NICMOS observation came 
from HST program 11118, L. Sromovsky PI, and was 
taken at 4:39 UT on 28 July 2007, at a phase angle of 
2.00° and a sub-observer latitude of 0.61°. The Keck 
II image was taken at 14:39 UT on 31 July 2007, at 
a phase angle of 1.87° and a sub-observer latitude of 
0.51°.] That these two observations sense roughly the 
same level in the atmosphere is demonstrated by the 
penetration depth plot in Fig. [6l which also displays the 
filter transmission functions. Although the Keck image 
has the higher spatial resolution, both images are ad¬ 
equate to resolve the gross latitudinal variations of in¬ 
terest. The I/F profiles for these two filters near the 
2007 Uranus equinox are displayed in Fig. 0 for /i=0.4 
and /i=0.6. These are plotted as true I/F values (not 
scaled in any way). At high latitudes in both hemi¬ 
spheres, and at both zenith angle cosines, the two pro¬ 
files agree with each other quite closely and are both 
increasing towards the equator. But at low latitudes, 
the reflectivity profile for the hydrogen-dominated wave¬ 
length continues to increase, while the profile for the 
methane-dominated wavelength decreases substantially, 
indicating that methane absorption is much higher at low 
latitudes than at high latitudes.This suggests that up¬ 
per tropospheric methane depletion (relative to low lat¬ 
itudes) was present at both northern and southern high 
latitudes in 2007, at least rough l y simi lar to the pattern 
that was inferred by iTice et ahl (j20I3[) from analysis of 
2009 IRTF SpeX observations. Latitudinal variations in 
aerosol scattering could distort these results somewhat, 
but because they affect both wavelengths to similar de¬ 
grees, most of the effect is likely due to methane varia¬ 
tions. 

4.2. Direct comparison of key STIS wavelength scans 

A similar spectral comparison of the 2012 STIS ob¬ 
servations can also be informative. By selecting wave¬ 
lengths that at one latitude provide similar I/F values 
but very different contributions by H 2 CIA and methane. 
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2007 Observations of Uranus 




Fig. 7.— Latitudinal profiles at fixed zenith angle cosines of 0.4 (A) and 0.6 (B) for F108N and PaBeta (Keck/NIRC2) filters taken near 
the Uranus equinox in 2007. At this point the southern hemisphere was still generally brighter than the northern hemisphere and the 38° S 
- 58° S southern bright band was still better defined and considerably brighter than the corresponding northern bright band. The relatively 
low equatorial I/F values for the methane-dominated PaBeta filter (1290 nm) indicates increased CH 4 /H 2 absorption at low latitudes. 


one can then make comparisons at other latitudes to see 
how I/F values at the two wavelengths vary with lati¬ 
tude. If aerosols did not vary at all with latitude, then 
this would be a clear measure of the ratio of CH 4 to 
H 2 . Fig. [ 8 ] displays a detailed view of I/F spectral re¬ 
gion where hydrogen CIA exceeds methane absorption 
(see Fig. [2] for penetration depths). Near 930 nm and 
827 nm the I/F values are similar but the former is dom¬ 
inated by methane absorption and the latter by hydro¬ 
gen absorption. Near 835 nm there is a relative mini¬ 
mum in hydrogen absorption, while methane absorption 
is still strong. At 50° N latitude and /i=0.6, I/F values 
are nearly the same at all three wavelengths, suggest¬ 
ing that they all produce roughly the same attenuation 
of the vertically distributed aerosol scattering. At low 
latitudes, as shown in Fig.[9l4, the I/F for the hydrogen- 
dominated wavelength increases, while the I/F for the 
methane-dominated wavelength decreases substantially, 
indicating an increase in the amount of methane rela¬ 
tive to hydrogen at low latitudes. Similar effects are 
seen for the 2002 observations. For yn=0.8 (Fig. [9^), 
which probes more deeply, the changes are even more 
dramatic. A color composite of these wavelengths (using 
R=930 nm, G= 834.6 nm, and B= 826.8 nm) is shown 
in Fig. [TOl where the three components are balanced to 
produce comparable dynamic ranges for each wavelength. 
This results in nearly blue low latitudes where absorption 
at the two methane dominated wavelengths is relatively 
high and green/orange polar regions as a result of the 
decreased absorption by methane there. 

The spectral comparisons in Figs. [9] and [10] also re¬ 
veal substantial secular changes between 2002 and 2012 . 
At wavelengths for which methane and/or hydrogen ab¬ 
sorption are important, the northern low-latitudes have 
brightened substantially, while the southern low latitudes 
have darkened. The 2002-2012 differences shown in Fig. 
[9] are not due to a view-angle effect because comparisons 
in that figure are made at the same view and illumina- 



Fig. 8 .— I/F and absorption spectra comparing the equilibrium 
H 2 CIA coefficient spectrum (divided by 1.2e-7, shown as dot- 
dash curve) and methane absorption coefficient spectrum (dashed). 
Note that the I/F spectrum has nearly equal I/F values at 826.8 
nm (A), 834.6 nm (B), and 930 nm (C), but H 2 absorption is 
much greater at A than at B, while the opposite is true of methane 
absorption, and at C only methane absorption is present. In a 
reflecting layer model, changes in cloud reflectivity should affect 
wavelengths A-C by the same factor, but changes in methane mix¬ 
ing ratio would affect C most and A least. 

tion angles for both years. The bright bands between 
38° and 58° (north and south) have also changed signif¬ 
icantly, with the southern band darkening dramatically, 
while the northern band brightened by a smaller amount. 
The northern bright band in 2012 was of lower contrast 
than the southern bright band in 2002. However, at con¬ 
tinuum wavelengths (see Fig. 4 in the analysis supple¬ 
ment file) secular changes are not evident. Aside from 
what appears to be a small calibration disagreement be¬ 
tween 2002 and 2012, in which corrected 2002 I/F values 
are about 2 % higher (at continuum wavelengths), the 
latitudinal variations are very similar for both epochs. 
[The corrected 2002 calibration we use here is based on 
WFPC2 comparisons and leads to I/F values that are 
3% smaller than published by KT2009.] The small size 
of these continuum differences is partly a result of the 
relatively smaller impact of particulates at short wave¬ 
lengths where Rayleigh scattering is more significant. At 
absorbing wavelengths the optical depth and vertical dis- 
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Fig. 9.— I/F vs latitude at /i = 0.6 (A) and fi = 0.8 (B) for three wavelengths with different amounts of methane and hydrogen absorption. 
Thin curves are from 2002 and thick curves from 2012. These are plots of center-to-limb fitted values instead of raw image data. In both 
cases the methane-dominated wavelengths have much reduced I/F at low latitudes. The 2002 calibration is based on WFPC2 comparisons 
and leads to I/F values 3% lower than the original KT2009 calibration. 


tribution of particulates have a greater fractional effect 
on I/F and thus small secular changes in these parame¬ 
ters can be more easily noticed. 

5. A SIMPLIFIED MODEL OF METHANE/HYDROGEN 
VARIATIONS 

Here we model detailed relative latitudinal variations 
in effective methane mixing ratios using a simplified fast 
spectral model that we put on an absolute scale with 
selected full radiative transfer model results. 

5.1. Model description 

KT2009 estimated the latitude variation of the CH 4 
volume mixing ratio using a simple model to fit the spec¬ 
tral region where hydrogen and methane have compara¬ 
ble effects on the observed I/F spectrum. They assumed 
that in the region from 819 nm to 835 nm the I/F spec¬ 
trum could be approximated by 

I/F = exp(Co + Cl X kcH^ + C 2 x kH 2 )^ (2) 

where kcH 4 and kH 2 are absorption coefficients of 
methane and hydrogen respectively. While this looks like 
a reflecting-layer model, it is actually more general. The 
coefficients Ci and C 2 give the weighted path lengths 
of observed photons. Because the methane coefficient 
is essentially independent of temperature and pressure, 
the weighting for coefficient Ci is the methane mixing 
ratio. But because kH 2 not only depends on temper¬ 
ature, but also on the square of the density, C 2 has a 
pressure-dependent weighting. Mean path lengths are 
also affected by aerosols. The observed variation of C 2 
is about 10 % peak-to-peak, aside from features that are 
consistent with noise. We expect that Ci would have a 
similar variation if the methane mixing ratio is constant 
with latitude, and if the path lengths of photons do not 
vary much across the selected spectral region. Because 
Cl and C 2 have different weightings, one might expect 
a little different variation for Ci, perhaps 5 % or 15 %. 
Thus, if the ratio C 1 /C 2 shows latitudinal variations of 
about 5 %, that could be due to aerosols. If they are 


much larger, it strongly suggests that they are due to 
a variation of the methane mixing ratio. Whether the 
path lengths vary across the spectral region or not can be 
tested by fitting the spectral shape with Eq. O The fits 
shown in Fig. [TT]are good. In fact, the slight misfit at the 
peak of the curves is likely due to a slight change of path 
lengths, indicating that the chosen spectral interval can¬ 
not be expanded without introducing systematic errors. 
To evaluate the hydrogen absorption coefficient, KT2009 
used an effective fixed temperature of 80 K, which we also 
used because it provided an approximate overall best fit 
to the spectra over a range of latitudes. We also followed 
KT2009 in assuming equilibrium hydrogen. 

While KT2009 chose ten 1 . 6 -nm wide segments from 
819 to 835 nm, we here extended the range down to 815 
nm to better sample the para hydrogen absorption peak 
near that wavelength, giving us 50 wavelengths over this 
range. The disadvantage of the wider interval is that fit 
quality is lower and somewhat more variable. However, 
sample calculations show very little difference in the de¬ 
rived latitudinal profile. In computing the estimates 
of fit quality, we assumed a somewhat arbitrary value of 
2% in relative I/F calibration-absorption modeling error 
and combined that with the center-to-limb fitting errors 
using the root of the summed squares. We used an ar¬ 
bitrary scaling of the kH 2 coefficients, then carried out 
a non-linear regression to obtain best fit model coeffi¬ 
cients, uncertainties, and the ratio of C 1 /C 2 and its un¬ 
certainty. Additional uncertainty is present due to the 
simplifications of the assumed model physics. The sam¬ 
ple fits provided in Fig. [TTlshow that the model is quite 
successful at fitting the observed spectra, and also allows 
tight constraints on the model coefficients and on the ra¬ 
tio of coefficients, which we use as a proxy for the CH 4 
/H 2 ratio. Section 5 of the on-line analysis supplemen¬ 
tal file provides a more extended discussion of this model 
and our attempts to constrain both effective temperature 
and hydrogen para fraction as well as relative absorber 
amounts. 
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Fig. 10.— Color composites of fitted center-to-limb smoothed images for 2002 (A, B, D) and 2012 (B-D) using color assignments R = 
930 nm (all methane), G = 834.6 nm (methane and hydrogen), and B = 826.8 nm (mostly hydrogen). North is up in all images except 
the reflected 2002 image in B. The blue tint at low latitudes for both years is due to locally increased methane absorption. In B the 2002 
image is inverted after remapping to the same central latitude as the 2012 image. The obvious asymmetry is not surprising because Uranus 
was in southern fall in 2002 and northern spring in 2012. In C we compare 2012 STIS observations with a Keck/N1RC2 near-lR image 
taken in the same year, high-pass Altered to enhance cloud structure. Note the small cloud features north of the high-latitude bright band. 
In D we compare 2002 and 2012 STIS observations at the same latitudes (i.e. without reflection of 2002 about the equator) and with both 
observations remapped to place the equator at the center of the image. In D the enhancement is the same for both 2002 and 2012, but 
using the 2012 albedo calibration function for both data sets (the solid curve in Fig. [d]). The southern high latitudes in 2002 were brighter 
and whiter than corresponding northern latitudes in 2012, providing evidence for a seasonal lag. The dark band near the center of the two 
images is 5° south of the equator. A comparison of latitudinal variations at fixed view angles is provided in Fig. [9] 




Fig. 11.— The 2012 I/F spectrum (solid curves) from 815 to 835 nm at ii=0.6 for two sample latitudes: 0° (left) and 65° (right), using 
the simplified model (+) of Eq. [2] as discussed in the text. Also shown are relative variations of absorption by methane (dotted) and H 2 
(dot-dash for para and double-dot-dash for ortho). 


5.2. Model results for latitude dependenee of CH 4 /H 2 . 

Fitting this crude model to every latitude for both 2002 
and 2012 center-to-limb fits leads to the latitude depen¬ 
dence shown in Fig. [121 where latitudes shown for each 
data set include the measured ranges and slight extrapo¬ 
lations of about 7° of latitude (the modest uncertainties 
of such projections are illustrated in Fig. 2 of the on-line 
supplement). The estimated error bounds shown here, 
which are twice as large as the formal fitting errors, are 
based on splitting the spectral data into two indepen¬ 
dent subsets and comparing the two profile results. Here 
we use a scaling of the CxjC^ ratio that best matches 
the methane volume mixing ratio (VMR) values esti¬ 
mated from full radiation transfer modeling that prop¬ 
erly accounts for the density dependence and temper¬ 
ature dependence of hydrogen absorption (discussed in 


subsequent sections). It is noteworthy that a single scale 
factor applied to all latitudes yields consistency with the 
full radiation transfer results at three different latitudes 
with different aerosol reflectivities. Using the same scale 
factor on the 2002 profile leads to about a 4% overall 
average between 30° S and 20° N. There are substantial 
differences between 2002 and 2012 where the two profiles 
overlap, though many of the small scale variations occur 
at nearly the same latitudes. The relative variations be¬ 
tween 35° S and 35° N have a correlation of 70%, with 
the 2012 variations being about 10% larger than the 2002 
variation. Yet the average 2012 methane VMR is closer 
to 4.5% at low southern latitudes where the 2002 profile 
is closer to 3.8%. Some of this difference might be due 
to the difference in aerosol structure, which is indicated 
by the higher reflectivity of the southern hemisphere in 
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Fig. 12 .— A: Latitude dependence of the ratio C 1 /C 2 , which is 
crudely proportional to the CH 4 /H 2 mixing ratio, scaled to best 
match the 2012 CH 4 VMR estimates at 10° S, the equator, and 
60° N, obtained from radiation transfer modeling. Those results 
are plotted as larger filled circles with error bars. The same scal¬ 
ing used on the 2002 data (thin dashed line) leads to a slightly 
lower methane VMR at low latitudes. Error bars for the ( 7 i /(72 
ratio points are provided only for scattered samples for clarity. B: 
The aerosol refiectivity term (exp(CO)), for 2002 and 2012, using 
the KT2009 spectral calibration divided by 1.03 for 2002 results. 
Dotted lines interpolate across the regions where bright bands are 
seen (38° - 58°). 

2002 (see Fig. [9j). Most of the small scale variability 
seen in the 2002 ratio is due to the hydrogen absorption 
term. Because this occurs without any evident small 
scale variability in the aerosol term (Fig.[9p), it suggests 
the possibility that small-scale temperature variations or 
vertical mixing variations that produce changes in para 
fractions might be responsible (this is discussed further 
in the next paragraph). A significant part of the small 
scale structure seen in the 2012 methane VMR values 
is also due to hydrogen coefficient variations, especially 
north of 60° N where the model cloud reflectivity varies 
relatively smoothly with latitude (Fig. [T^i. at least at 
these wavelengths (small variations are seen at near-IR 
wavelengths, as discussed later). 

Also noteworthy is the difference between 2012 north 
and 2002 south polar regions (there is no overlap between 
2002 and 2012 in either polar region). In 2002 the mixing 
ratio from 75° S to 50° S is close to 2.4% and exhibits 
only small variations with latitude, while the 2012 pro¬ 
file at equivalent northern latitudes displays substantial 
and nearly sinusoidal variations with a latitudinal wave¬ 
length near 10°. It should be noted that each of these 
points is derived from 50 wavelengths and is thus far less 
noisy than one might suspect based on the plots of single 
wavelengths, as in Figs. [9] and ini The reality of these 
high-latitude features is further supported by the fact 
that the same analysis applied to similar latitudes in the 
opposite hemisphere using 2002 observations did not find 
such large variations, suggesting that the variations are 
more likely due to a hemispheric difference rather than an 
inherent high latitude artifact of the analysis. Also, there 


is no obvious reason why such artifacts should suddenly 
disappear at latitudes below 60° N. Further evidence sup¬ 
porting the reality of these latitudinal variations in the 
methane to hydrogen absorption ratio is the discovery 
of a correlated pattern of cloud reflectivity variations in 
high-resolution Keck imagery, which is discussed in the 
analysis supplement file. A hint of this pattern can also 
be seen in Fig.fTOlO. We found similar patterns in August 
and November images and both had moderately signifi¬ 
cant negative correlations with the methane VMR (^ 9 % 
for August and(^ 5-10 % for November). While the cor¬ 
relations are suggestive, it is still possible that the STIS 
latitudinal variation is due to noise and new observations 
will be needed to confirm this feature. The explanation 
for this feature, presuming that it is real, might be that 
latitude bands of reduced above-cloud methane lead to 
clouds appearing locally brighter due to reduced methane 
absorption. An alternate explanation is that the local 
variation in methane to hydrogen absorption is due to 
changes in the para fraction of hydrogen. If we force the 
methane mixing ratio to be constant over the 55° - 82° 
latitude range and fit spectral variations by adjusting the 
para fraction and effective temperature, then para frac¬ 
tion increases appear where we previously had methane 
decreases. The para increases suggest downwelling and 
thus we have a correlation of descending motion with re¬ 
duced cloud amount (or reduced cloud reflectivity). This 
alternative interpretation presents a more appealing pic¬ 
ture than our initial interpretation, although it has a 
slightly worse spectral fit quality. 

6 . FULL RADIATIVE TRANSFER MODELING OF 
METHANE AND AEROSOL DISTRIBUTIONS 

The spectral profile comparisons (Figs. [71 [HI and [9]), 
and especially the somewhat more quantitative model¬ 
ing of Fig. [T 2 I strongly suggest that there is a permanent 
and more or less symmetric high latitude methane deple¬ 
tion on Uranus. But, because hydrogen absorption has a 
density squared dependence, and the methane has a ver¬ 
tically varying mixing ratio, a more accurate constraint 
on methane requires full radiative transfer modeling, in¬ 
cluding the effects of more realistic aerosol distributions, 
as described in the following. The more complete model¬ 
ing is also needed to provide the scaling factor that con¬ 
verts the C 1 /C 2 ratio plotted in Fig. [12] to a methane 
VMR. 

6.1. Radiation transfer calculations 

We used the r adiation transfer code described by 
iSromovs^ (|2QQ5aD , which include Raman scattering and 
polarization effects on outgoing intensity, though this is 
a minor virtue at the wavelengths employed in our anal¬ 
ysis. To save computational time we emplo yed the ac¬ 
curate polarization correction described bv ISromovskvl 
(l2nn5]nH . After trial calculations to determine the effect 
of different quadrature schemes on the computed spec¬ 
tra, we selected 14 zenith angle quadrature points per 
hemisphere and 14 azimuth angles. Calculations with 
10 quadrature points and 10 azimuth angles changed fit 
parameters by only about 1%, which is much less than 
their uncertainties. To characterize methane absorption 
at CCD wavelengths we used the coefficients of KT2009. 
To model collision-induced opacity of H 2 -H 2 and He-H 2 
interactions, we interpolated tables of absorption coef- 
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ficients as a function of pressure and temperature that 
were com puted with a prog:ra m provided by Alexandra 
Borysow (jBorvsow et al.ll2QQQf ). and available at the At¬ 
mospheres Node of NASA’S Planetary Data System. We 
assumed equilibrium hydrog:en for mo s t calc ulations, fol¬ 
lowing KT2009 and iSromovskv et al.l (j2Qll[ ). 

6.2. Cloud model 

The cloud model we used (referred to as the compact 
model in what follows) is a modification of the KT2009 
diffuse model (containing 4 vertically extended layers) 
that is fully described in Section 7 of the supplemental 
analysis file. Our five-layer model is illustrated in Fig. 
M and summarized in Table O In this table, parameter 
names begin with m for layers containing Mie particles 
and with hg for layers containing particles with double 
Henyey-Greenstein phase functions. Here we use Mie 
particle scattering as a convenient parameterization of 
wavelength dependence, realizing that the particles are 
generally not spherical. Adjustable parameters are iden¬ 
tified in the Value column of the table. We chose this 
model because it provided substantially better fits than 
the diffuse model. Our top two haze layers use the same 
parameterization as KT2009. The main change we made 
to the KT2009 model was to replace its middle tropo¬ 
spheric layer with two compact layers: an upper middle 
tropospheric cloud layer (UMTC) and a lower middle 
tropospheric cloud layer (LMTC). The UMTC layer is 
composed of Mie particles, which we characterized by a 
gamma size distribution with an adjustable mean parti¬ 
cle radius and a fixed normalized variance of 0.1, a fixed 
refractive index of 1.4, and an imaginary index of zero. 
For the LMTC layer we used particles with the same 
scattering properties as given by KT2009 for their tropo¬ 
spheric particles. Both of these compact layers have the 
bottom pressure as a free (adjustable) parameter and a 
top pressure that is a fixed fraction of 0.93 x the bottom 
pressure. This degree of confinement is approximately 
the same as obtained for the c loud layer inferred from 
the radio occultation analysis (jSromovskv et al.ll201l[ ). 
The motivation for introducing these replacement layers 
was to obtain more flexibility in vertical structure and 
to allow the possibility of including a thin cloud near the 
methane condensation level, as suggested by the occul¬ 
tation analysis. 

The last change we made was to replace the KT2009 
bottom tropospheric layer by a compact cloud layer at 
5 bars (the BTC, or bottom tropospheric cloud), with 
adjustable optical depth and with the KT2QQ9 tropo¬ 
spheric scattering properties. ISromovskv et al.l (j2Qlll ) 
found that this layer was needed to provide accurate fits 
near 0.56 and 0.59 /im, but its pressure was not well con¬ 
strained by the observations (pressure changes could be 
compensated by optical depth changes, to produce es¬ 
sentially the same fit quality). Whether this deep cloud 
is vertically diffuse or compact also could not be well 
constrained. Sample phase functions and the scatter¬ 
ing properties of particles used in our cloud model are 
plotted vs. wavelength in the analysis supplement file 
(Section 7). Note that the KT2009 empirical functions 
defining wavelength-dependent single-scattering albedo 
and phase functions only apply over the 300 nm to 1000 
nm range for which they were derived. The role played 
by each model layer in creating a spectrum that matches 


0.01 


CO 0.1 

cc 

-Q 

(D 

3 


10 

0.001 0.01 0.1 1 10 100 
Optical depth per bar 

Fig. 13.— Compact cloud model vertical structure. Compared to 
the KT2009 model, this has two layers with additional adjustable 
parameters of base pressure (for the new hg2 layer) and base pres¬ 
sure and particle radius (for the new m2 layer). 

sample observations is also provided in the analysis sup¬ 
plement file (Section 8). 


ml 


. hg1 

Compact Cloud Model 

m2, m2_p, m2_r CH4 cloud; 

- 

hg2, hg2_p *^28? 

hgS NH4SH?- 



6.3. Fitting procedures 

To avoid the complexity of fitting a wavelength- 
dependent imaginary index in the methane layer (the 
UMTC layer in Table [2]) we fit only the wavelength range 
from 0.55 /im to 1.0 /im. We chose a wavenumber step 
of 118.86 cm“^ for sampling the observed and calculated 
spectrum. This yielded 69 spectral samples, each at three 
different zenith angle cosines (0.3, 0.4, and 0.6), for a 
total of 207 points of comparison. Our compact layer 
model has seven to eight adjustable parameters (see Ta¬ 
ble [2]), leaving 201 or 200 degrees of freedom {Np) re¬ 
spectively. (The expected value of is equal to Np for 
an optimized fit, and the formal uncertainty in any fit 
parameter is given by the change in its value that in¬ 
creases by one unit when the remaining parameters 
are optimized (jPress et al.lll99^ .) We fixed the BTC 
base pressure to 5 bars (see Table [2]) and usually ignored 
the UTH layer because of its negligible effects. We used 
a modi fied Levenberg-Marquard t non-linear fitting algo¬ 
rithm (|Sromovskv and Frvll201Cih to adjust the fitted pa¬ 
rameters to minimize and to estimate uncertainties in 
the fitted parameters. Evaluation of y^ requires an esti¬ 
mate of the expected difference between a model and the 
observations due to the uncertainties in both. We fol¬ 
lowed t he same approach for est imating uncertainties as 
used by ISromovskv et al.l (|2Qllf ) , which combined mea¬ 
surement noise (estimated from comparison of individual 
measurements with smoothed values), modeling errors of 
1%, relative calibration errors of 1% (larger absolute cal¬ 
ibration errors were treated as scale factors), and effects 
of methane absorption coefficient errors, taken to be ran¬ 
dom with RMS value of 2% plus an offset uncertainty of 
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TABLE 2 

Compact cloud model parameters 



Parameter / function 

Value 

Stratospheric haze 
of Mie particles 
with gamma size 

mljpd (bottom pressure) 
ml_r (particle radius) 
mlJ) (variance) 

fixed at 100 mb 
fixed at 0.1 pm 
fixed at 0.3 

distribution 

nl (refractive index) 
ml^odpb (optical depth/bar) 

nr=1.4, ni given by KT2009 
adjustable! 

Upper tropospheric haze 
of double HG particles 
(UTH) 

hgl-pd (bottom pressure) 

'Uj(X) (single-scatt. albedo) 
phase function (KT2009) 
hgl^odpb (optical depth/bar) 

fixed at 0.9 bars (top pressure = ml-pd) 
KT2009 function 

gi = 0.7, g2 = 0.3, /i(A) given by KT2009 
adjustable2 

Upper middle tropospheric 
cloud layer of Mie 
particles (UMTC) 

m2_p (bottom pressure) 
m2_r (particle radius) 
m2_b (variance) 
n2 (refractive index) 
m2 od (optical depth) 

adjustable3 (top pressure = m2_px0.93) 

adjustabled 

fixed at 0.3 

fixed at n=1.4 

adjustableS 

Lower middle tropospheric 
cloud of double HG 
particles (LMTG) 

hg2jp (bottom pressure) 
vj{X) (single-scatt. albedo) 
phase function (KT2009) 
hg2 od (optical depth) 

adJustableO (Ptop = hg2jp X 0.9) 
given by KT2009 

gi = 0.7, g2 = 0.3, /i(A) given by KT2009 
adjustable7 

Bottom tropospheric cloud 
(BTC) 

ro(A) 

phase function (DHG) 
hgS-od (optical depth) 
hg3jp (bottom pressure) 

same as previous layer 
same as previous layer 
adjustableS 
fixed at 5 bars 


NOTE: KT2009 equations defining wavelength dependent parameters are reproduced in the analysis sup¬ 
plement file. Usually hgl-odpb was found to be too small to bother including in our fits. 


5x10“^ (km-amagat)“^. The uncertainty in is ex¬ 
pected to be ^25, and thus fit differences within this 
range are not of significantly different quality. 

7. COMPACT MODEL FITS 

The compact model described previously provides gen¬ 
erally better fits than the diffuse models, reaching 
values in the low to mid 200’s compared to mid 300’s 
to low 500’s for the diffuse model. This is presumably 
due to its greater flexibility to fit vertical distributions 
because the main cloud layer is divided into two sublay¬ 
ers with adjustable pressures as well as adjustable optical 
depths. Another difference from diffuse models is that its 
upper sublayer has its wavelength dependence controlled 
by Mie scattering, which is parameterized by particle ra¬ 
dius and refractive index, rather than the wavelength- 
dependent double Henyey-Greenstein parameters of of 
the lower sublayer, which we take to be those of KT2009 
and reproduced in the analysis supplement file. 

In the following, we first constrain the effective 
methane mixing ratio profile at key latitudes to define the 
scaling factor for Fig. [121 fhen try to fit cloud structure 
to the observed spectra over a wide range of latitudes, 
using the best-fit equatorial methane vertical profile to 
discover where it does not apply. Poleward of 30° N we 
obtain improved fits by reducing the effective mixing ra¬ 
tio. But the best fits are obtained with a depletion profile 
that exhibits decreasing depletion with increasing depth, 
which also is more physically plausible than a uniform 
depletion at all depths. In other words, the preferred 
latitudinal profile of methane becomes more uniform at 
greater depths. 

7.1. Constraining equivalent methane mixing ratios at 
key latitudes 

The detailed latitudinal variation of the equivalent 
methane mixing ratios plotted in Fig. [12] is based on a 


scaled ratio of fit coefficients [CxjC^. To define that 
scaling, we used fits that fully account for the differ¬ 
ent vertical absorption profiles of hydrogen CIA and 
methane, carried out at three latitudes covering the range 
of mixing ratio variation (10° S, 0°, and 60° N). To esti¬ 
mate the optimum equivalent deep methane mixing ratio 
at each latitude we did compact model fits for a variety 
of methane profiles (Dl, DE, El, EE, El, EG, and G) 
that have a range of deep methane mixing ratios (2.22%, 
2.76%, 3.20%, 3.6%, 4.00%, 4.5%, and 4.88%, respec¬ 
tively). These profiles are all cons istent with the Voy - 
ager 2 occultation measurements of iLindal et al.l (|1987[) . 
and thus also have different temperatur e and above-cloud 
metha ne profiles. They are all from I Sromovsky et al.l 
(|2011[ ) except for DE and EG, which we constructed to 
be occultation consiste nt using the same techniques as 
I Sromovsky et al.l (|2011[ ). The fit results at our selected 
key latitudes are plotted in Eig. [TH where we show two 
fit quality measures as a function of the deep methane 
mixing ratio of the profiles that were employed. The two 
measures are (1) and (2) the signed error at 825 nm 
divided by the error expected just from measurement er¬ 
rors. At 10° S the minimum x^ occurs at a lower methane 
mixing ratio (3.6%) than the minimum error at 825 nm 
(4.6%), while at the equator the two measures agree on 
4%. For 60° N, we find by extrapolating that the 825- 
nm error reaches zero at 1.6%, while the corresponding 
X^ minimum is between 2.2% and 2.6%. We consider the 
825-nm estimate to be the most relevant since it provides 
the most direct comparison between methane and hydro¬ 
gen absorptions and uses the /i = 0.6 spectrum to weight 
the deep atmosphere more than the overall fit, which uses 
/i = 0.3, 0.4, and 0.6. However, its greater uncertainty 
limits led us to use intermediate estimates between these 
two measures: 4.1±0.5% at 10° S, 4.0±0.5% at the equa¬ 
tor, and 2.0±0.5% at 60° N. These points were plotted 
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Fig. 14.— Best-fit compact model parameters vs deep methane 
mixing ratio for un depleted occultation-cons istent Dl, El, EF, FI, 
and G profiles from ISromovskv et al.l (|201in . and DE and FG pro¬ 
files created using the same procedure. The fits are to STIS spectra 
at 10° S, the equator, and 60° N, for zenith angle cosines of 0.3, 
0.4, and 0.6. The values are shown with solid curves, and the 
825-nm error ratio to expected error at /r=0.6 is shown with dotted 
curves. Horizontal arrows indicate which axis to read in each case. 

in Fig.[T 2 ]as filled circles. That these results lead to such 
a consistent scaling for the simplified model justifies its 
use in estimating latitudinal variations. The low-latitude 
average mixing ratio is sufficiently close to the FI profile 
value of 4% that we used FI as our base profile. 

Note that an occultation consistent profile with a deep 
mixing rati o as low a s 1.6% , which is the deep VMR 
assumed by iTice et aTI (j2Q13f ) , would be grossly inconsis- 
tent with the obs erved spectra at low latitudes. Though 
iTice et aTI (j2Q13f ) were able to fit their spectra with such 
low methane VMR values, they may have been able to 
make up for lacking methane absorption at longer wave¬ 
lengths by assuming a low single-scattering albedo in 
their main cloud layer. Also, based on sample spectral 
fits shown in their Figs. 4, 7, and 15, they did not fit very 
well the critical spectral region near 825 nm, where their 
model spectra show relatively more hydrogen absorption 
than the observations, which suggests that their methane 
mixing ratio is too low. They also reduce the weight of 
this spectral region in their fits by assigning relatively 
large uncertainties compared to other regions with simi¬ 
lar I/F values. Their fits were also aided by using much 
more above-cloud methane than would be allowed by the 
0% relative h umidity implied by the Voyager 2 occulta¬ 
tion results of iLindal et al.l (j 19871 ) for their assumed deep 
mixing ratio. 

7.2. Compact model fits versus latitude for the FI 
profile 

We first tried to fit the compact model to a range of lat¬ 
itudes, while assuming that the FI profile of temperature 
and methane applied to all latitudes. The entire set of 
results is plotted in Fig. [151 results for latitudes up to 
30° N are presented in Table O There are several remark¬ 
able features of these fits. First, the ml vertically diffuse 
stratospheric haze layer is found to have significantly re¬ 
duced opacity north of 25° N. Whether this is a real 
effect and a real change from the structure inferred from 
2002 observations or a result of a different spectral cali¬ 
bration using WFC3 observations is unknown. Perhaps 


a related difference is that the upper sublayer (UMTC) 
of the main cloud layer, a Mie layer, has a fitted particle 
radius generally in the 0 . 2 - 0 .5 micron range, compared to 
the 1 . 2 -micron radius th at seemed to be preferred by the 
iSromovskv et al.l (j 2011 f ) analysis of 2002 STIS observa¬ 
tions. This is likely due to the revised albedo calibration 
function (Fig. |4]), which produces a stronger I/F decline 
with wavelength that is better matched by smaller par¬ 
ticles (see analysis supplement file). 

At the equator the results are also unusual. In diffuse 
model fits we found that the upper tropospheric haze 
(the second layer in Table [2]) provided a sharply increased 
contribution at the equator, and seemed to be the main 
change associated with the bright equatorial band visible 
in many images at wavelengths of intermediate methane 
absorption (e.g. Fig. |3K, and S). In the compact model 
fits, however, changes in the next lower and higher layers 
(first and third layers in Table [2]) seem to be relatively 
more important, with changes in pressure and particle 
size of the third layer (m 2 _p and m 2 _r) as well as the 
optical depth per bar of the top layer {ml-odpb) playing 
a role. The derived value of the optical depth per bar 
of the second layer (hgl-odpb) for the compact model is 
comparable to its uncertainty and was ignored for most 
compact model fits. 

Between 30° S and 30° N, the 825-nm signed error, 
given by (model-observed)/uncertainty, is relatively fiat. 
The small positive error seen in this at most latitudes in 
this range suggests that the mixing ratio in the FI pro¬ 
file is not quite large enough. The rather large negative 
825-nm errors found at high latitudes are associated with 
excessively high methane mixing ratios in the model cal¬ 
culations. The measure of fit quality is generally low, 
but is a more erratic function of latitude over this range. 
Beyond 40° N, the 825-nm signed error becomes much 
larger as the overall fit quality measured by becomes 
much worse. This indicates that the assumed latitude- 
independent profile of methane is not consistent with 
the observations. A different methane profile is clearly 
needed at high latitudes as was already demonstrated at 
60° N from fits used to define the scaling of the simplified 
model. 

7.3. Why methane depletion should not extend to great 

depths 

While structure and methane profiles with a 2 % CH 4 
deep mixing ratio provide the best fits at 60° N com¬ 
pared to other profiles with constant mixing ratios up to 
the condensation level, the deep contrast between equa¬ 
tor and pole (2% VMR to 4% VMR) is not physically 
plausible. Such deep latitudinal gradients in composition 
would lead to density gradients along isobars. As a con¬ 
sequence of geostrophic and hydrostatic bala nce, these 
gradi ents would lead to vertical wind shears (jSun et al.l 
EMU). These vertical wind shears acting over the entire 
atmosphere would likely lead to cloud top winds that 
would be highly incompatible with the observed winds. 
Thus we considered methane vertical distributions in 
which the higher latitudes had depressed mixing ratios 
restricted to shallow depths in the upper troposphere 
only. As indicated by KT2009, the 2002 spectral observa¬ 
tions did not requ ire that methane d epleti ons extend to 
great depths, and ISromovskv et ah] (j 2011 f ) showed that 
shallow depletions were preferred by the 2002 spectra. 
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TABLE 3 

Best-fit parameters for equatorial compact cloud layer models versus latitude assuming simulated FI profile for 
LATITUDES BETWEEN 30° S AND 30° N, AND USING DEPLETED FI PROFILES FOR HIGHER LATITUDES. 


Lat. 

o 

m2_p 

bar 

hg2jp 

bar 

ml-odpb 

m2_od 

hg2_od 

hg3-od 

m2_r 

pm 


(m-o)/u 

-30 

1.19T0.03 

1.41T0.02 

0.21T0.04 

0.34T0.09 

0.96T0.05 

4.09T0.75 

0.32T0.12 

299 

-0.32 

-20 

1.31T0.03 

1.51T0.04 

0.25T0.04 

0.46T0.11 

0.94T0.07 

3.64T0.90 

0.19T0.01 

373 

0.69 

-10 

1.25T0.02 

1.52T0.03 

0.42T0.04 

0.39T0.09 

1.10T0.07 

3.84T0.90 

0.18T0.03 

338 

0.85 

0 

1.18T0.02 

1.48T0.03 

0.62T0.03 

0.43T0.07 

1.19T0.06 

3.94T0.67 

0.47T0.05 

205 

-0.07 

10 

1.27T0.02 

1.47T0.03 

0.34T0.04 

0.47T0.10 

1.12T0.07 

4.83T1.09 

0.16T0.02 

359 

1.00 

20 

1.26T0.02 

1.47T0.03 

0.26T0.04 

0.48T0.10 

1.16T0.07 

5.42T1.20 

0.19T0.03 

344 

0.29 

30 

1.19T0.02 

1.36T0.02 

0.09T0.04 

0.24T0.04 

1.01T0.04 

5.05T0.72 

0.20T0.02 

227 

0.52 

38 

1.24T0.03 

1.58T0.03 

0.23T0.03 

0.37T0.07 

0.95T0.05 

2.84T0.44 

0.47T0.06 

204 

-0.23 

45 

1.23T0.04 

1.53T0.03 

0.21T0.04 

0.38T0.09 

1.06T0.06 

1.86T0.43 

0.45T0.08 

307 

-0.18 

50 

1.17T0.05 

1.57T0.02 

0.09T0.05 

0.27T0.04 

1.07T0.04 

1.07T0.24 

0.22T0.02 

226 

0.29 

60 

1.24T0.05 

1.70T0.04 

0.12T0.04 

0.36T0.06 

0.82T0.03 

0 

0.43T0.05 

234 

0.14 

70 

1.32T0.05 

1.85T0.07 

0.04T0.04 

J 

0.48T0.08 

0.65T0.04 

0 

0.43T0.05 

246 

0.33 


methane volume mixing ratio of 4%, and a He/H 2 ratio of 0.1306. At other latitudes the fits used methane 
depleted profiles with vx=2>.Q and depletion depths between the value that minimizes 825-nm errors and those 
that minimize (both listed in Table |4] and plotted in Fig. I15l> . The uncertainty in x^ is ~25 and thus fits 
differing by less than this are not of significantly different quality. The column labeled (m — 6)1 u is the fit 
error at 0.825 fim expressed as the ratio of (model I/F - observed I/F) to the estimated I/F uncertainty for fi 
= 0.6. These results are plotted in Fig. nni The parameter hgl^odpb was omitted from these fits for reasons 
discussed in the main text. 



Fig. 15.— Best-fit compact model parameters vs. latitude assuming the FI structure and methane mixing ratio profile (left) and using 
the best-fit depletion depth profiles (right). Error bars indicate fitting uncertainty. In each case, the lower panel displays fit quality in 
terms of x^ (+) and signed 825-nm error (o). In the bottom left plot, the behavior of the 825-nm error indicates that the actual methane 
absorption declines below the assumed value beginning about 40° N. In the bottom right plot, both 825-nm error and x^ benefit from use 
of depleted profiles. In the bottom panels horizontal arrows indicate which axis to read for each curve. 
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We will show here that the 2012 spectral constraints favor 
relatively shallow methane depletions at middle latitudes 
but increasingly deeper depletions at high latitudes. 

Here we constructed profiles with shallow CH 4 deple- 
tion using the “proport ionally descended gas” profiles of 
(jSromovskv et al"]l 2 Qll[ ) in which the model FI mixing 
ratio profile a{P) is dropped down to increased pressure 
levels P'{ol) using the equation 

P' = P X [1 + {a{P)/adr{Pd/Pc - 1)], (3) 

for Ptr < P < Pd^ where Pd is the pressure depth at 
which the revised mixing ratio o:'{P) = ol{P') equals 
the uniform deep mixing ratio ad^ Pc is the methane 
condensation pressure before methane depletion, Ptr is 
the tropopause pressure (100 mb), and the exponent vx 
controls the shape of the profile between 100 mb and Pd- 
The profile s with vx = 1 are similar in fo rm to those 
adopted by iKarkoschka and Tomaskol (j2011f ) . 

7.4. Compact model fits at 6(F N vs depletion depth 

To account for the lower methane mixing ratio that 
is observed at high latitudes we need to find a methane 
sink that can create local depletions. The most obvi¬ 
ous one is methane condensation, which can cause the 
low-temperature above-cloud region to have much lower 
methane mixing ratios than the regions below. When 
this depleted gas is mixed downward, it can cause lo¬ 
cal depressions of the methane mixing ratio in the upper 
troposphere. Thus we expect only the upper troposphere 
to be depleted and a significant question is whether we 
can constrain the depth of that depression using STIS 
observations. 

Here we start with an FI profile at 60° N and then use 
Eqn. [3] to deplete the methane above a base pressure Pd 
at a rate controlled by a secondary parameter vx, which 
causes rapid depletion with decreasing pressure when set 
to small valu es and smaller depletion rates when set to 
larger values. iSromovskv et al.l (j 2 Qll[ ) found that vx = 1 
seemed to be preferred by the 2002 data set, but we found 
larger values preferred by the 2012 data set (this might 
also be true for 2002 if explored more extensively and 
with the 2012 calibration function). We tried a range 
of Pd values for vx ranging from 1 to 5. Subsets of our 
results are plotted in Fig. [161 which displays and 825- 
nm error values as a function of Pd for different assumed 
values of vx. The optimum values are summarized in 
Table jH We found that minima appear at increasing 
depths as the depletions become more gradual with de¬ 
creasing pressure (for larger values of vx). The best fit at 
60° N was found for ux = 3.0, with Pd = 30 bars to min¬ 
imize 16 bars to minimize the 825-nm error. For 

vx = 2, the minimum is near 10 bars and somewhat 
larger than the minimum for vx = 3, while the minimum 
825-nm error is found at 7.5 bars, closer to the mini¬ 
mum than for the vx = 3 case. Clearly, ux = 1 is a very 
poor choice. But even that profile fits better (y^ = 320) 
than the D1 profile (y^ = 360), which is the best fitting 
of the occultation consistent profiles with vertically uni¬ 
form mixing ratios. The vertical variation of methane 
VMR for the best fit Pd value is shown for each vx case 
in Fig. [TTK. Note that shallower and deeper depletions 
all produce similar mixing ratios near 1.7 bars, and that 
greater depletions at depth result in somewhat higher 



Fig. 16.— Fit quality estimates for compact cloud layer model 
fits to spectra at 60° N (x^ on left axis, 825-nm error on right), 
shown as function of depletion depth for depletion profiles with 
different vx values. Depletions are relative to an FI base model. 

mixing ratios at pressures below 1.7 bars. It should also 
be noted that the STIS measurements aren’t sensitive to 
the large depths implied by the best-fit Pd values. It is 
fitting the shape of the empirical profile at much lower 
pressures that constrains the Pd value. Thus it would 
be worthwhile to explore other depletion functions with 
different vertical variations. 

7.5. Latitudinal variation in depletion depth and aerosol 

structure 

To fix the problem of poor spectral fits at latitudes 
beyond 30° N, we chose vx = 2 and vx = 3 vertical 
variation functions and found the best-fit value of Pd as 
a function of latitude for each. At all latitudes we found 
that vx = 3 provided the best fits. The y^ and 825-nm 
error results as a function of Pd for vx = 3 are plotted in 
Fig. 12 of the analysis supplement. Here we list best-fit 
Pd values in Table [Hand the best-fit aerosol parameters in 
Table |3l where the aerosol model is chosen for Pd values 
between those giving the smallest y^ and those giving the 
smallest 825-nm error. The best-fit aerosol results are 
also plotted in Fig. [15] (right panels) for comparison with 
the fits using the undepleted FI profile at all latitudes. 
The depletion depth increases from 30° N to 70° N (the 
highest latitude fit) and in this region y^ and 825-nm 
errors have been greatly reduced relative to the prior fits 
with undepleted profiles. 

One puzzling result is that the cloud layer that we 
thought was associated with methane condensation (the 
UMTC or m2 layer in Fig. continues into the polar 
regions where methane condensation should not occur 
because of the decreased methane mixing ratio. Perhaps 
we should not have used the same labels for m2 points 
at high latitudes that we used at latitudes from 30° S 
to 30° N (see right panel of Fig. [15]). In fact, the sharp 
drop in optical depth of the m2 layer between 20° N and 
30° N might have continued northward, and what we 
identified as the m2 layer north of 30°, might actually 
be a continuation of the hg2 cloud. In that case, the 
layer we identified as the hg2 cloud in these northern 
fits might actually have a different composition from the 
corresponding layer at low latitudes. Also note that the 
hg3 layer declines dramatically north of 50° N, where Pd 
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Fig. 17. — Depletion models for 60° N with different vx values 
vx=3.0 for all latitudes. 

exceeds the assumed depth of the cloud layer. The polar 
region possibly has a completely new set of cloud layers. 

8. DISCUSSION: METHANE AS A CONSTRAINT ON 
MERIDIONAL MOTIONS. 

Fig. [TSl illustrates a mechanism for depleting methane 
at high latitudes that is essentially the same as what 
was previously suggested by KT2009. Rising methane- 
rich gas at low latitudes is dried out by condensation and 
sedimentation of methane ice particles. That dried gas is 
then transported to high latitudes, where it begins to de¬ 
scend, bringing down methane depleted gas, which then 
gets mixed with methane-rich gas on its return flow. The 
figure illustrates two different styles of return flow. With¬ 
out lateral eddy mixing across the streamlines, the only 
restoration of the methane mixing ratio would be through 
evaporation of the precipitating condensed methane at 
low latitudes. The depth of the depletion at high lat¬ 
itudes might be controlled by the depth of the merid¬ 
ional cell or the depth at which cross-streamline mixing 
predominates. The suggested circulation would promote 
formation of optically thin methane clouds (or hazes) 
at low latitudes, but inhibit methane cloud formation at 
high latitudes. This seems consistent with the lack of ob¬ 
served discrete cloud features south of 45° S. However, 
we hav e seen discrete cloud fea tures at high northern lat¬ 
itudes (|Sromovskv et al]l2QQ9[) , which might very well be 
composed of something other than methane, given that 
their cloud tops se em to be deeper than the methane 
condensation level (jSromovskv et alll 2 Q 12 [ ). The general 
downwelling would also tend to inhibit all condensation 
clouds, as sub-condensation mixing ratios would be cre¬ 
ated by such motions. However, this does not rule out 
localized regions of upwelling and formation of conden¬ 
sation clouds occupying a small fractional area. 

The circulation cell in Fig. [18] would need to be very 
deep to be consistent with microwa ve observations prob¬ 
ing the 10-100 bar region of U ranus (|de Pater et al .11 19891 
119911: iHofstadter et al.l I2QQ7I ) . These reveal a symmetry 


and best fit depletion models for different latitudes (B), assuming 


S Pole 45 S equator 45 N N Pole 

methane depleted by 
condensation 



Fig. 18.— Speculative meridional flows in which rising flows at 
low latitudes produce methane condensation (white region) that 
depletes the gas of methane; the depleted gas is transported to high 
latitudes where it descends to reduce the mixing ratio at higher 
pressure. The depletion depth might be limited by the depth of 
the circulation, or by lateral mixing. A very deep cell of this type 
has been inferred from microwave observations (see main text). 
Note that the streamlines of the northern cell provide an example 
of shallower depletions at intermediate high latitudes than those of 
the southern cell, the former option being more consistent with the 
observed increasing methane depletion depth at higher latitudes. 

pattern in which microwave absorbers (NH 3 , H 2 S) are 
depleted at both high southern and high northern lati¬ 
tudes, suggesting a non-seasonal equator-to-pole merid¬ 
ional circulation, with upwelling at low latitudes and 
down-welling at high lat i tudes (Ide Pater and Lissanerl 
IMTntlde Pater et, al JlT^ ITMil lHi^t,a,rlt,er et, a,lJl2fin7r 

similar to the circulation cells illustrated in Fig. [T 8 | If 
this deep cell extended to the upper troposphere, it could 
be consistent with a shallow CH 4 depletion, as long 
as the flow at the several bar level was dominated by 
poleward flow that did not go through the drying-by¬ 
condensation process, as suggested by the inner stream¬ 
line in the northern flow pattern in Fig. [T 8 | While our 
best fits indicate that the largest fractional depletions 
occur in the upper troposphere, at higher latitudes some 
degree of depletion could extend deeper than the 10 -bar 
level. Because these deep meridional cells would likely be 
dominated by deep atmospheric conditions, they would 
probably have the same symmetry properties as the deep 
atmosphere (symmetry about the equator), which would 
suggest qualitatively that the north and south polar re- 
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TABLE 4 

Optimum depletion depth based on fit quality and 825-nm error as a function of depletion profile rate (small vx provides 

SHARPER DEPLETIONS ABOVE THE DEPLETION DEPTH LIMIT Pd)- 


P lanet ogr aphic 
Latitude, ° N 

vx 

Vd 

minimum 

825-nm error 

minimum 825-nm error 

Pd, bars 825-nm error 

38 

2.0 

2.0 

215 

-0.50 

3.0 

0.28 

226 

38 

3.0 

3.0 

204 

-0.23 

4.0 

O.IO 

207 

45 

2.0 

2.0 

318 

-0.77 

3.0 

-.07 

326 

45 

3.0 

4.0 

307 

-0.18 

4.0 

-0.18 

307 

50 

2.0 

4 

239 

-O.OI 

4 

-0.01 

239 

50 

3.0 

10 

225 

0.57 

8 

0.29 

226 

60 

1.5 

5.5 

282 

0.19 

5.0 

0.00 

282 

60 

2.0 

10 

238 

0.48 

7.0 

-0.12 

252 

60 

3.0 

30 

227 

0.60 

16 

0.01 

242 

60 

5.0 

200 

243 

0.34 

90 

-0.02 

249 

70 

2.0 

12 

253 

0.50 

8 

-0.18 

263 

70 

3.0 

30 

246 

0.33 

20 

-0.09 

254 


gions should not look very different. With respect to 
the gross character of the methane depletion, they don’t 
look very different, even though they do exhibit different 
cloud morphologies. 

A problem with a single deep cell in each hemisphere 
is that it would seem to inhibit formation of H 2 S con¬ 
densation clouds that are good candidates for producing 
the bright bands that form between 38° and 58° in both 
hemispheres. In fact, such a simple circulation would 
tend to inhibit all condensation clouds at high latitudes, 
in spite of the fact that cloud particles of some kind are 
detected there. A 3-layer circulation cell (Fig. [19]) can 
provide more opportunities for condensation cloud for¬ 
mation. In this case the deep cell produces ammonia de¬ 
pletion through formation of an NH 4 SH cloud at low lat¬ 
itudes, an intermediate cell produces high latitude con¬ 
densation of H 2 S, and the top cell provides high-latitude 
depletion of methane. This structure is also qualita¬ 
tively compatible with the observed equatorward motion 
of the long-lived and largest d iscrete cloud feature seen 
on Uranus, known as the Berg (iSromoyskv and Frvll2QQ5l: 
iSromovskv et al.l I2QQ9I: Ide Pater et al.ll 2 Qll[) . This fea¬ 
ture is thought to have been at the level of the putative 
H 2 S cloud, and in this model would ride along with the 
equatorial ffow near the 1.5 bar level. However, the re¬ 
semblance between the speculated structure in Fig. [T9l 
and the structure shown in the right panels of Fig. [15] is 
crude at best. While there seems to be a change in cloud 
structure between high and low latitudes, if the interpre¬ 
tation is restricted to models with condensation clouds 
only, it is hard to explain the existence of what seems to 
be an H 2 S cloud at low latitudes where the three-layer 
cell structure would inhibit such a cloud. Furthermore, 
the drift of the Berg may not be a relevant constraint 
on meridional ffow. If the cloud features comprising the 
Berg were generated by an unseen vortex, its drift may 
be controlled more by the vorticity of the zonal flow th an 
by weak meridional flows (|Lebeau and Dowlindll99§) . 

An additional problem with the above suggested 
meridional circulation structures is that they do not ap¬ 
pear to be consistent with Voyager observations. Using 
thermal emission measurements by th e Voyager IRIS in¬ 
strument in 1986 iFlasar et al.l (|I987[) inferred upwelling 
motions at ±30°, where relatively lower temperatures 


S Pole 45 S equator 45 N N Pole 



Fig. 19.— Speculated 3-layer cell structure providing equator- 
ward flow at the berg level and poleward flow above methane con¬ 
densation and NH 4 SH cloud formation levels. Here cloud forma¬ 
tion/condensation regions are indicated by light rectangles. 

were measured, and subsidence at both equator and 
poles. The temperatures were measured in the 60-200 
mbar and 600-1000 mbar regions, and the circulations 
inferred from a zonally symmetric linear model with fric¬ 
tional and radiative damping. But this symmetric pat¬ 
tern is no t evident in th e para fraction distribution in¬ 
ferred by iConrath et al.l (jI998[) . which indicates a sea¬ 
sonal contrast between hemispheres at the time of the 
Voyager encounter, with the northern hemisphere having 
lower temperatures and higher para fractions, consistent 
with radiative cooling and downward motion. Hadley cell 
configurations with hemispherical symmetry are also not 
consistent with recent n umerical modeling of t he seasonal 
circulation of Uranus bv iSussman et all (| 2 QI 2 [) . which in¬ 
dicates that there should be tropospheric cross-equatorial 
flow peaking near equinox. Another point worth noting 
is that the presumed upwelling at low latitudes, which is 
supported by the inference of an optically thin methane 
cloud in those regions, is not accompanied by an abun¬ 
dance of discrete cloud features at those latitudes. Like¬ 
wise, the observed cloud structure at high northern lati¬ 
tudes, which appear somewhat deeper than expected for 
condensed methane, may be more related to localized 
upwellings rather than any broadly defined meridional 
cells. This may be similar to the situation on Jupiter 
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in which a region of overall downwelling (a belt) con¬ 
tains many examples of localized c onvection and light¬ 
ning (jShowman and de Paterl[2QQ5f ). It is also conceiv¬ 
able that local increases in methane humidity might af¬ 
fect aerosol opacity even though the methane mixing ra¬ 
tio remained below saturation. The effect of water vapor 
humidity on the size and scattering properties of hygro¬ 
scopic aeroso l is significant and well known in the ear th’s 
atmosphere (jPilat and CharlsonlfTQGGl: lKastenlll969l ). It 
is conceivable that a haze of hydrocarbon aerosols orig¬ 
inally formed in the stratosphere might have their scat¬ 
tering properties altered in regions of varying methane 
humidity, even without the occurrence of methane con¬ 
densation. Such an effect might even give the appear¬ 
ance of a condensation cloud in regions of local upwelling. 
However, we do not know if the background aerosols on 
Uranus would interact with methane in the same way 
as hygroscopic aerosols interact with water vapor in the 
earth’s atmosphere. 

Neptune provides an interesting second example of 
methane depletion at high latitudes. Using 200 3 STIS 
spectra of Neptune, iKarkoschka and Tomaskol (j2QIIf ) 
showed that between 80° S and 20° N the variation in 
Neptune’s effective methane mixing ratio was very sim¬ 
ilar to that observed for Uranus. To explain this dis¬ 
tribution they also suggested a mechanism similar to 
what is illustrated in Fig. [181 In this case, the 45° - 
50° S latitude band, one of t he most active reg ions of 
cloud formation on Neptune (jKarkoschkal 120111 ) turns 
out to be a region where downwelling is needed to ex¬ 
plain the methane depletion results. However, recent pa¬ 
pers suggest a meridional circulation for Neptune that 
has upwelling in this region. From a re-analysis of Voy- 
ager/IRIS 25-50 /im mapping of tropos pheric temper¬ 
atures and para-hydrogen disequilibrium iFletcher et al.l 
(j2014l ) suggested a symmetric meridional circulation with 
cold air rising at mid-latitudes and warm air sinking at 
equator and poles. Based on a multi wavelength anal¬ 
ysis t hat included near-IR to microwave observations in 
2003. Ide Pater et HI (j2014j ) also detected warm polar and 
equatorial regions, where they infer downwelling motion, 
and cooler middle latitudes, where they infer upwelling 
motion. Such a circulation pattern is inferred to extend 
to great depths and would seem to be in conflict with 
the pattern needed to produce the observed upper tro¬ 
pospheric methane depletion. Thus, both planets seem 
to have more complicated stories than we are currently 
able to explain with simplified models. 

9. SUMMARY AND CONCLUSIONS 

We observed Uranus with the HST/STIS instrument 
in 2012, aligning the instrument’s slit parallel to the spin 
axis of Uranus and stepping the slit across the face of 
Uranus from the limb to the center of the planet, build¬ 
ing up half an image with each of 1800 wavelengths 
from 300.4 to 1020 nm. The main purpose was to con¬ 
strain the distribution of methane in the atmosphere of 
Uranus, taking advantage of the wavelength region near 
825 nm where where hydrogen absorption competes with 
methane absorption and displays a clear spectral signa¬ 
ture. Our analysis of STIS observations of Uranus from 
2012 and comparisons with similar 2002 observations, as 
well as analysis of imaging observations from 2007, have 
led us to the following conclusions. 


1. In 2012 the methane mixing ratio in the upper tro¬ 
posphere of Uranus was depleted at high northern 
latitudes (relative to equatorial values), especially 
beyond 40° latitude, to a degree very similar to 
what was inferred for southern high latitudes from 
2002 STIS observations. This is based on (1) di¬ 
rect spectral comparisons of STIS latitudinal pro¬ 
files at wavelengths with similar penetration depths 
but different amounts of hydrogen absorption, (2) 
simplified quantitative modeling of the 815-835 nm 
spectrum, and (3) full radiative transfer modeling. 

2. We also found that the north and south depletions 
wer e simultaneous in 2007, as also suggested by 
the iTice et al ] (l2nT^ analysis of 2009 SpeX cen¬ 
tral meridian spectra, and thus probably not a sea¬ 
sonal effect. Our 2007 result is based on direct 
spectral comparisons near the equinox, using an 
HST/NICMOS F108N image that is sensitive to H 2 
absorption and a Keck/NIRC2 PaBeta-filtered im¬ 
age that senses about the same atmospheric level, 
but is dominated by methane absorption. 

3. We followed KT2009 in using a simplified model of 
the 815 nm - 835 nm spectral region to estimate 
the relative latitudinal variation of the methane 
volume mixing ratio at 1° intervals. When this 
relative variation was absolutely scaled to match 
effective mixing ratios determined by full radiative 
transfer modeling at 10° S, the equator, and 60° 
N, we found that the effective mixing ratio varied 
on both large and small spatial scales. At the large 
scale we found the VMR to increase from roughly 
2% within 30-40° of the poles to about 4% within 
20-30° of the equator. The 2012 observations sug¬ 
gest an overall increase in methane VMR at low 
latitudes by about 0.5% relative to 2002. However, 
this might be due to a change in the vertical dis¬ 
tribution of aerosols at low latitudes rather than a 
change in methane VMR. 

4. Between 60° N and 82° N, the simplified model 
revealed a variation of ±0.005 in methane VMR 
relative to a local mean of 0.020, and a nearly si¬ 
nusoidal variation with a period of about 10° in 
latitude. A similar variation was not seen in south 
polar regions using the same analysis techniques on 
a similar STIS data set from 2002. Nevertheless, 
there is a chance that this feature is not real and 
needs confirmation by further measurements. 

5. Using Keck2/NIRC2 high-pass filtered H-band 
images from 16 August and 4 November 2012, 
we computed zonally averaged brightnesses and 
compared their variations with latitude with the 
methane VMR variations between 55° N and 82° 
N. We found similar patterns in the latitudinal vari¬ 
ations of cloud reflectivity in August and Novem¬ 
ber images and both had significant negative cor¬ 
relations with the methane VMR modulations in 
the same latitude region, suggesting that latitude 
bands of reduced above-cloud methane make the 
clouds appear brighter due to reduced methane 
absorption. An alternative interpretation is that 
the methane mixing ratio is relatively constant in 
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this region, with para fraction variations (associ¬ 
ated with local vertical circulation cells) explain¬ 
ing the variation in relative strengths of CH 4 and 
H 2 absorptions. In this case downwelling correlates 
with reduced cloud reflectivity. 

6 . At 60° N, we tried a variety of vertical variation 
functions for the methane depletion profile, charac¬ 
terized by the exponent ux, which we varied from 
1 to 5. We found that the depletion depth Pd in¬ 
creased as the sharpness of the depletion decreased, 
with best-fit mixing ratios near 1.7 bars being com¬ 
parable for all choices of vx. The best fit of all was 
obtained for vx = 3, with vx = 2 also providing a 
good fit with somewhat better agreement between 
the Pd that minimized (10 bars) and the Pd that 
minimized the 825-nm error (7.5 bars). 

7. We carried out fits to determine depletion depth as 
a function of latitude, assuming a depletion profile 
shape defined hy vx = 2 and also by ux = 3, with 
the latter shape providing the best fit at all lati¬ 
tudes. We found that Pd increased with latitude, 
beginning just beyond 30° N. At high latitudes the 
depth of the downwelling flow could exceed ten bars 
or more, although we are only sensitive to the upper 
tropospheric depletions (P < 10 bars, as indicated 
by the penetration depth plot of Fig. 2), so that 
the relatively large depletion depths we found may 
be partly a result of the particular empirical func¬ 
tion we used in our models. Other profile shapes 
might be able to fit the data without producing as 
great a depth of depletion. 

8 . Using the depleted profiles to constrain the aerosol 
parameters, we found a lowering in the altitude (in¬ 
crease in base pressure) of the m2 layer north of 30° 
N (the third aerosol layer from the top). Since the 
lowered methane mixing ratio also implies that this 
layer can no longer be associated with widespread 
methane condensation (at lower latitudes it is lo¬ 
cated at the methane condensation level), it might 
here be composed of other materials or produced 
by widely dispersed local convective events, or pro¬ 
duced by changes in background aerosols due to 
absorption of methane instead of condensation of 
methane. 

9. The association of high-latitude methane deple¬ 
tions with descending motions of an equator-to- 
pole deep Hadley cell does not seem to be consistent 
with the behavior of the detected aerosol layers, at 
least if one ignores other cloud generation mech¬ 
anisms such as sparse local convection. Both on 
Uranus and Neptune, aerosol layers seem to form 
in what are thought to be downwelling regions on 
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the basis of the effective methane mixing ratio de¬ 
terminations. A three-layer set of circulation cells 
offers some advantages in producing condensation 
clouds, but also fails to provide a good match to 
the detected aerosol layers. 

10. Secular changes between 2002 and 2012 vary with 
wavelength. At continuum wavelengths changes 
appear to be very small. When the KT2009 cali¬ 
bration is adjusted by 3% to match WFPC 2 band¬ 
pass filtered images, the 2002 I/F value are found 
to be 2 % greater than corresponding 2012 values. 
This residual difference may be due to calibration 
uncertainties. On the other hand, at wavelengths 
with noticeable gas absorption (as in Fig. [9]), the 
northern hemisphere has brightened considerably 
since 2002 , by about 25% at mid latitudes at 827 
nm, and the southern hemisphere has darkened, by 
about 25% at mid latitudes at 827 nm. 

In the future, better constraints on the vertical profile 
of methane as a function of latitude could be addressed 
by additional modeling work with the 2012 STIS spec¬ 
tra, trying different functional forms for vertical deple¬ 
tion profiles. More detailed analysis at more latitudes 
using full radiation transfer modeling for both 2002 and 
2012 would be useful in clarifying whether low-latitude 
changes between 2002 and 2012 are real. Additional 
STIS observations in future cycles are also needed to con¬ 
firm the northern high-latitude modulations we found in 
the apparent mixing ratio. Additional quantitative con¬ 
straints might also be derived from analysis of the verti¬ 
cal wind shears that are implied by the horizontal density 
gradients associated with latitudinal compositional gra¬ 
dients. Additional work with numerical circulation mod¬ 
eling might also be productive in understanding how the 
methane mixing ratio affects and is affected by atmo¬ 
spheric circulation patterns. 
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SUPPLEMENTAL MATERIAL. 

An expanded discussion of analysis techniques and results is provided in the file ura_stis_analysis_supplenient.pdf, 
which can be found on line at http://www.ssec.wisc.edu/planetary/uranus/onlinedata/ura2012stis/, The hy- 
perspectral cube containing calibrated I/F values as a function of wavelength and location, with a naviga¬ 
tion back plane that provides viewing geometry and latitude-longitude coordinates for each pixel, is also pro¬ 
vided there in ura_stis2012_datacube.fits. A detailed explanation of the contents of this file can be found in 
the file README_SUPPLEMENTAL.TXT. A sample IDL program that reads the cube file, plots a monochro¬ 
matic image, extracts data from a particular location on the disc, and plots a spectrum, is provided in the file 
stis_cube_read_example.pro. The IDL astronomy library will be needed to read the data cube, which is in the EITS 
format. 


